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We present an effective operator formalism for open quantum systems. Employing perturbation theory and 
adiabatic elimination of excited states for a weakly driven system, we derive an effective master equation which 
reduces the evolution to the ground-state dynamics. The effective evolution involves a single effective Hamil- 
tonian and one effective Lindblad operator for each naturally occurring decay process. Simple expressions are 
derived for the effective operators which can be directly applied to reach effective equations of motion for the 
ground states. We compare our method with the hitherto existing concepts for effective interactions and present 
physical examples for the application of our formalism, including dissipative state preparation by engineered 
decay processes. 
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I. INTRODUCTION 

Understanding the dynamics of open quantum systems |[T) 
requires methods to model their temporal dynamics. Gener- 
ally the evolution of an open quantum system is described by 
the master equation which determines the evolution of the sys- 
tem's density matrix p. Solving the master equation for the 
full density matrix is in many situations cumbersome. This is 
also the case for closed systems where the evolution can be 
described by wave functions. Since in general an open system 
involves both unitary and dissipative dynamics, the complex- 
ity of its evolution is, however, substantially higher than for 
the corresponding closed system, as the number of entries in 
the density matrix is the square of the number of entries in the 
wave function. To make the description of an open quantum 
system manageable as well as to gain physical insight into the 
evolution of the system it is therefore desirable to develop ef- 
fective theories which reduce the complexity of the system. 

In an open system there are often quite different time scales 
associated with different effects, such that the Hilbert space 
can be divided into two parts, one for the rapidly decaying 
(excited) states, and one for the comparably stable (ground) 
states. For instance, for weakly driven atoms the evolution and 
decay of the excited states happen on a time scale which is fast 
compared to any other time scale in the system. In such situ- 
ations it is desirable to eliminate the rapidly evolving excited 
states to get a simpler description of the slow evolution of the 
ground states. A standard method for doing this is adiabatic 
elimination [2], where the density matrix equations involving 
the excited states are solved by assuming a slow evolution of 
the ground states. This can then be used to describe an ef- 
fective evolution of the ground states. This procedure can, 
however, be rather involved, as there are many density matrix 
elements if the system is large. 

In this paper we present a simple method to eliminate the 
excited states and to reduce the system dynamics to the ground 
states. The method we present is essentially equivalent to adi- 
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abatic elimination but is much easier to apply in practice. By 
formalizing the procedures leading to adiabatic elimination 
we obtain simple expressions for the effective operators de- 
scribing the ground-state evolution. With our expressions one 
avoids the often tedious steps leading to adiabatic elimination 
and can obtain the effective operators by evaluating simple 
formulas. In particular, we have found that these methods are 
very convenient for studying dissipative state preparation (3J- 
13, where the goal is to engineer decay processes such that a 
system evolves into a desired state. To this end it is highly 
desirable to have a convenient tool to rapidly identify the ef- 
fective dissipative dynamics of the system. 

For closed systems with purely unitary couplings similar 
simplifications have been used to derive effective Hamiltoni- 
ans, in particular by James and co-workers ll8i iT0l . If decoher- 
ence is added to the system, the joint unitary and dissipative 
dynamics can be captured by introducing non-Hermitian de- 
scriptions. The use of complex energies allows for the combi- 
nation of the energy of a resonance with its width [11]. Corre- 
spondingly, non-Hermitian Hamiltonians are commonly used 
to describe the dynamics of open systems |[T2"14T6l . In quan- 
tum optics the use of non-Hermitian Hamiltonians is put in a 
more rigorous form by the so-called quantum jump formal- 
ism, or Monte Carlo wave function method lfT2l . which is 
equivalent to the evolution by a Markovian master equation. 
In this method the non-Hermitian Hamiltonian describes the 
evolution in the absence of decay, whereas quantum jumps are 
introduced at random times to account for the resulting state 
after a decay. 

The separation of the Hilbert space into rapidly and slowly 
evolving ground and excited states, similar to what we con- 
sider here, has been studied also for non-Hermitian Hamilto- 
nians. In particular, for a coupling of the ground to the excited 
states much weaker than the evolution inside the subspaces, 
a formalism for effective processes is provided by the Fesh- 
bach projection-operator approach JT3- Similar formalisms 
for successive nuclear reactions have been studied by Wei- 
denmuller and co-workers ff8l . Methods based on the Fes- 
hbach projection operator method have been used in several 
fields [19|. However, these treatments are only concerned 
with the evolution from the effective non-Hermitian Hamil- 
tonian and therefore ignore the quantum jumps describing the 
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FIG. 1: (Color online) Ground and excited subspaces and couplings. 
The nondecaying ground states (corresponding to projector P g ) are 
coupled to the decaying excited states (projector P c ) by the perturba- 
tive (de-) excitations V+ (VI) [solid lines connecting the subspaces 
(red)]. The Lindblad operators Lk represent various decay process 
(dotted lines). The couplings inside the ground and excited subspaces 
are given by Hg [solid lines inside the lower subspace (green)] and 
H c [solid lines inside the upper subspace (blue)], respectively. 



state after the decay. The procedure we present here is an ex- 
tension of the Feshbach projection-operator approach to also 
include these quantum jumps. As a result, our formalism can 
be used to describe the full evolution of the density matrix of 
the system after elimination of the excited states. This gener- 
alization is crucial for describing situations where we are also 
interested in the state of the system after a decay. 



II. AIM OF THIS PAPER 

Before proceeding we first outline the main result of this 
paper. As illustrated in Fig. [T] we assume the open system 
to consist of two distinct subspaces, one for the ground states 
and one for the decaying excited states. The couplings of these 
two subspaces are assumed to be perturbative. Furthermore, 
we assume that the dynamics of the system is Markovian such 
that the time evolution of the density operator p can be de- 
scribed by a master equation of Lindblad form 
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with effective Hamilton and Lindblad operators 
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connecting only the ground states. Here V+ (V-) are the per- 
turbative (de-) excitations of the system and H g is the ground- 
state Hamiltonian. H^h 1S the non-Hermitian Hamiltonian of 
the quantum jump formalism 
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with H c being the Hamiltonian in the excited-state manifold. 
The effective master equation of Eq. Q provides an approx- 
imation of the dynamics in Eq. ([T]i by the effective dynamics 
of its ground states. Thus, the effective operator formalism 
allows for a substantial reduction of the complexity of the dy- 
namics of an open system. In essence, adiabatic elimination of 
the excited states in the presence of both coherent and dissipa- 
tive processes is formalized in a compact form by the effective 
operators of Eqs. ([3]l and (HI. 

In the following section, Sec. [Hi] we will derive the effec- 
tive operators of Eqs. Q and Q. In Sec. HIE we show a first 
elementary application of our formalism to a driven dissipa- 
tive two-level system. Thereupon, we discuss the possibility 
to use our effective operator formalism to engineer decay pro- 
cesses in a four-level system in Sec. Ill F 

In Sec. IIV Al we turn to an extension of the formalism 
which includes nonperturbative ground-state couplings. This 
extended formalism will then be used for a detailed analysis 
of the effective processes in a three-level Raman system in 
Sec. IIV Bl In Sec. IIV CI we introduce a second extension of 
the effective operator formalism that allows for several per- 
turbations or fields. The most general formalism is presented 
in Sec. 
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A comparison to similar existing methods is 
provided in Sec. [V] Readers who are more interested in the 
applications than in the derivations may turn to Sees 
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III F and IV B where we give examples of effective coher- 
ent and dissipative processes in typical quantum systems and 
discuss a simple dissipative state preparation scheme. 



III. DERIVATION OF THE EFFECTIVE OPERATORS 



where H is the Hamiltonian of the system and each of the 
Lindblad operators L k represents a source of decay which we 
assume to take the system from the excited to the ground sub- 
space. By combining perturbation theory of the density oper- 
ator and adiabatic elimination of the excited states we reduce 
the dynamics to an effective master equation involving only 
the ground-state manifold 



We now present the derivation of the effective equation of 
motion in Eqs. Q-Q. The evolution of the density opera- 
tor pit) in the Schrodinger picture is governed by the master 
equation of Lindblad form, given by Eq. ([TJ. The Hamiltonian 
H stands for unitary couplings of the system, such as coher- 
ent driving. Each Lindblad "jump" operator L% accounts for 
a dissipative process, such as spontaneous emission. 
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A. Projection-operator formalism 

We use the projection-operator method of Feshbach ifTTl 
to structure the Hilbert space into two subspaces, one for the 
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ground states and one for the excited states, represented by 
the projection operators P g and P c , with P g + P c = 1 and 
P g P c = 0. Accordingly, we divide the Hamiltonian into four 
parts: 

H = H g + H c + V++V-. (6) 

Here, the interactions inside the ground subspace are labeled 
by H g = PgHP g , and inside the excited subspace by H c = 
P C HP C . The perturbative excitations V+ = P c HP g and de- 
excitations V- = PgHPc (V+ = V- and V = V+ + V-) 
connect the two subspaces. 

We assume the ground states as stable and the excited states 
to be decaying to the ground states. The Lindblad operators 
can then always be written as = PgL^P a . The mentioned 
couplings inside and between the subspaces are illustrated in 
Fig.0 

B. Non-Hermitian time evolution in the quantum jump picture 

Combining unitary and dissipative dynamics within a single 
non-Hermitian Hamiltonian has widely been studied in vari- 
ous areas of physics, as mentioned in the Introduction. In 
quantum optics the use of non-Hermitian Hamiltonians is for- 
malized by the so-called quantum jump picture |[T2l . in which 
an effective non-Hermitian Hamiltonian describes the evolu- 
tion of the system in the absence of a quantum jump. In order 
to distinguish the non-Hermitian Hamiltonian of the excited 
states from the effective (Hermitian) Hamiltonian of Eq. (|3]l 
we have denoted it as i?NH in Eq. (jij. It incorporates the 
excited state Hamiltonian H c and the decay terms of the anti- 
commutator part of the master equation (flj. Introducing -Hnh 
to Eq. ([T]i we obtain a reduced master equation 

P = - i ((£nh + H g + V)p - p(fr NH + H e + V )) 
+ Y,L k pL\ (7) 

k 

Here, we have included the decay terms which describe the 
loss of population from the excited states to 7?nh in the com- 
mutatorlike Hamiltonian part. The last "feeding term" which 
describes the gain of the population of the ground states by 
decay from the excited states remains. 

For ground-state interactions H g much weaker than those 
between the excited states H c , the dynamics of the decay- 
ing excited states are mainly governed by the non-Hermitian 
Hamiltonian 7?nh- As all excited states are decaying, all 
eigenvalues of i?NH are nonzero so that its inverse exists 
within the excited-state subspace. 

C. Perturbation theory in the interaction picture 

In the following we assume the couplings of the ground and 
excited subspaces V± to be sufficiently weak to be described 



as perturbations of the evolution governed by an unperturbed 
Hamiltonian H n = H g + H^n- Based on this assumption we 
perform perturbation theory of the density operator. To this 
end we change into the interaction picture by a transformation 
with the operator 

6(f) = e~ i6ot = e -KHnn+B s )t _ (g) 
Then the reduced master equation of Eq. |7]i transforms into 

p{t) = - i (v(t)p(t) - p(t)vHt)) + LkWWtit) 

k 



(9) 

with operators transformed accordingly 

p(t) = d-\t)p(d- 1 {t))\ (io) 

V(t) = 0-\i)(H Q + V)0(t) + i^-O(t) 

= 6- 1 (t)vd(t), (ii) 

L k (t) = 6- 1 {t)L k 6{t). (12) 



To derive the effective operators we perform a perturbative 
expansion of the density operator in a small parameter e 

P(t) - ^(P i0) (t) + eP (1) (i) + e 2 /5 (2) (t) + ...) (13) 

and obtain a recursive formulation of the reduced master equa- 
tion in powers of e, 

p( n \t) = - ziV^p^it) - p^-VtyVi (t))+ (14) 

+ J2Lk(t)p (n) (t)Ll(t), (15) 

k 

where we have used that V is a small parameter V oc e. The 
first three orders of the recursive reduced master equation read 

^\t)=Y J L k {t)p W {t)L\{t), (16) 

k 

^Ht) = -i(v(t)pW(t)-p(V(t)vi(t)) + 

+ J2Lk(t)p W (t)Ll(t), (17) 

k 

P ) (t)=-i(y(t)p^{t)-p^{t)V\t)) + 

+ J2Lk(t)p^(t)Ll(t). (18) 

k 

In the absence of initial excitations, decay processes can be 
neglected for orders n < 1 so that 

P (o) W=0, (19) 

j>W(t) = - 1 (v(t)pW(t) - p^wHt)) » (20) 

p\t)=-i(y(t)p^{t)-p^{t)V\t)) + 

+ Y J Lk{t)p {2) {t)L\{t). (21) 

k 
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We use the projection operator approach for the density oper- 
ator to separate the evolution of ground and excited states. In 
doing so we reduce the evolution of the ground states to 

Pj^(t)P s =P e ^(t)P g = 0, (22) 
P gP ^(t)P g = - iP g (v(t)p^(t) ~ P (1) W\t)) P g + 
+ Y J Lk{t)Pj 2) {t)P c L\{t). (23) 

k 

In the last line we have used that for each Lindblad operator 
we can write = P g LkP c , as decay only occurs from the 
excited to the ground states. Consequently, the ground states 
are connected by unitary and dissipative processes of second 
order. Also note that since the transformation in Eq. ([8]) is 
nonunitary, the perturbation V(t) is non-Hermitian. For the 
dynamics of the excited states we find 



Pj^(t)P e = P e ^(t)P e = 



(24) 



P c p\t)P c = -iP c (v(t)pW(t) - pV(t)V*(tj) P c . 

(25) 

As we have assumed that the excited states do not gain pop- 
ulation from decay, Eq. (25 i does not exhibit any dissipa- 
tive feeding terms. Hence, the evolution of the excited uni- 
tary dynamics is solely driven by the interaction Hamiltonian 
V(t). While the dynamics of the second-order terms connect 



the states either in the ground or in the excited subspace we 
note that interactions between the subspaces are given by the 
first-order terms P g p (1) (t)P e and P e p (1) (t)P s - 



D. Adiabatic elimination of the excited states 

In principle, a solution to the remaining second-order mas 
ter equations for the ground and the excited states in Eqs. (|23 



and (25 i can be computed. This solution can, however, still be 



very complicated. In particular, if the decaying excited states 
are almost unpopulated, it is preferable to obtain a more com- 
prehensible solution. In the following, we choose to reduce 
the complexity of the dynamics by restricting it to the ground 
states. To this end, we perform adiabatic elimination of the 
excited states: 



P e j>W(t)P e nO. 



(26) 



Consequently, the dynamics of second order in Eq. (21 1 are 
approximated by the dynamics of the ground states given by 
Eq. Q23) . Below we follow the recursion of the perturbative 
expansion and carry out the perturbation inte gral s. 

We obtain P e p$(t)P e by integrating Eq. {25}, and 
by integrating Eq. 
into Eq 



(23 1 



20p , and insert the resulting expressions 
Having excluded the dynamics of the excited 
states by adiabatic elimination, we find the open system to 
evolve according to 



p g ^\t)p g = - p g v(t) Qf* dt> v-(f )p (( V)) p g - p g (7* dt' pW(t')vHt')) v\t)p g + 

+ P g Y J Lk{t)P c [ dt ' J* dt" {y{t^ Q \t ll )V\t ll ) + V{t l Tp^\t ll )V\t 1 )) P c L\(t)P g . (27) 



Here, we have omitted terms where the density operator is 
sandwiched between perturbations P g V and VP g . As p^> 
lives in the ground-state subspace, these terms do not con- 
tribute to the ground-state evolution and can therefore be ne- 
glected. The remaining expression in Eq. ( 27 1 contains two 



Hamiltonian-like and two Lindblad-like terms, for which we 
will carry out the integrals: 



h = P s V(t) / dt'V{t')~^\t')P g 



(28) 



h = P c [ dt' f dt" V(t')pW(t")V\t")P e . (29) 
Jo Jo 

In this section we assume the direct interactions within the 
ground-state subspace to be perturbative. Hence, the ground- 
state evolution is negligibly small compared to the one for the 
excited states so that we have 0(t)P g ~ P g . Consequently, 



Ii simplifies to 



h « V-6(t) (J dt' O- 1 ^')) V+p {0) (t) (30) 



Carrying out the integral we find 

(zi? NH ) _1 e^ NHt ' 



V + P i0) (t) 



V- 



(31) 



In the last step we have used that the term emerging from the 
lower limit of the integral at t' = maintains its time depen- 
dence of e" lff ™ ! and is therefore detuned with respect to the 
term originating from the integral limit at t' — t, i.e., by an 
approximation similar to the rotating wave approximation we 
keep the unity term in the expression 1 — exp(— iH-^nt). This 
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condition is equivalent to the standard approximation of adia- 
batic elimination and is justified provided that the time evolu- 
tion of the ground states is slow compared to the time scale set 
by ^?nh- The second term of Eq. (27 1 is treated accordingly, 
yielding the Hermitian conjugate of the result in Eq. pT) . 

For the last two Lindblad-type terms in Eq. (27 1 we carry 
out the double integral 12- To this end we approximate 



p (0) {if') in Eq. ( 29 1 by p^ {t). This can be argued the follow- 



ing way: Above we have assumed that the density matrix of 
the ground states p(°> (t) evolves slowly and to second-order 
in V. Another dependence on V{t') and V(t") would only 
involve features of fourth order in the evolution. We neglect 
these higher orders by dropping the dependence of p(°\t') 
and p^{t") on the time scales of V{t') and V{t"), which 
yields p^ {t). Thus, we can separate the integral and write 



1 



NH 



(32) 



Again we have assumed that the ground states are slowly vary- 
ing compared to the time scale of H^h so that 0(t)P g ~ P g 
and discarded detuned terms. The remaining term of Eq. (|27]> 



yields the same result as Eq. ( 32 1. 

We insert Eqs. ( 3 1 1 and ( 32 1 back into Eq. ( 27 1 and trans- 



form back into the Schrodinger picture. In doing so, we obtain 
the effective unitary and dissipative dynamics of the ground 
states, 

pj 2) p s = -i(A tB -i £( £ -) t£ k) p {0) + H - c -+ 

+ £i c V 0) (4 fe ff ) t . 03) 

k 

with an effective Hamiltonian and effective Lindblad opera- 
tors as defined in Eqs. ([3]l and Q. To reach this form we have 
used the equality 
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From here it can easily be seen that Eq. p3\ is equivalent to 
Eq. p|. Thus, we have reduced the unitary and dissipative 
dynamics of the open quantum system described by Eq. ([TJi to 
the effective master equation of Lindblad form in Eq. (|2]), ob- 
taining the effective Hamiltonian H e g and Lindblad operators 
L£ ff ofEqs. <U and 

The effective Hamiltonian of Eq. Q is the same as the orig- 
inal result of Feshbach ifTTl . In addition, we have found effec- 
tive Lindblad operators for second-order decay processes. As 
can be seen from Eq. (|4]), each of them consists of weak co- 
herent excitation V + , evolution between the excited states by a 
"propagator" ff^, and subsequent decay L^. Thus, adiabatic 
elimination of the excited states of an open quantum system is 
formalized in a compact manner by Eqs. d2T 
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FIG. 2: Two-level system. A ground state |0) is coherently coupled 
to an excited state 1 1) by a field with a Rabi frequency of f2. The 
excited state 1 1} is subject to spontaneous decay at a rate of 7. 



E. Example: The two-level system 

The most elementary configuration our formalism can be 
applied to is given by a two-level system with a ground state 
|0) coherently coupled to a spontaneously decaying excited 
state |1) (Fig. [2|. This example is rather straightforward 
also without the theory developed here; yet, we include it to 
demonstrate the general formalism. The Hamiltonian for this 
system reads 



H = H e + V 

H e = A\l)(l\,H g = 

K = §(|0}<1| + |1>(0|) 



(34) 
(35) 

(36) 



Here, |0) and 1 1) are coupled with a Rabi frequency of ft and a 
detuning of A. We write the (de-) excitation as V+ = ^ |1) (0| 
(V- = §|0)(1|). Spontaneous emission from the excited to 
the ground level at a rate of 7 is represented by the Lindblad 
operator, 

L = V7|0><1|. (37) 
Consequently, the non-Hermitian Hamiltonian is found to be 

&nh=A|1)<1|, (38) 

with a complex energy A = A — ^ that combines the detun- 
ing A and the decay width 7 of the excited state 1 1) to a single 
complex quantity A. 

By applying the effective Hamiltonian formula of Eq. ([3]) 
we adiabatically eliminate the excited state |1) and obtain an 
effective Hamiltonian for the ground state |0): 




+ H.c 
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"4A 2 +7 



, _ s |0><0| 
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(39) 



This effective Hamiltonian describes an effective ac Stark 
shift A c ft of level |0) caused by the coherent driving. By ap- 
plying Eq. Q together with V, 7?nh, and L 7 as specified 
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FIG. 3: (Color online) Dissipative state preparation in a four-level 
system. An effective decay process consisting of weak excitation 
from ground state \gi) with a resonant Rabi frequency of £7, evolu- 
tion between the excited levels |ei) and [62), coupled with g, and 
subsequent decay k prepares ground state \g%). The effective decay 
rate K e g is engineered by the choice of the detunings of the excited 
states, A and 8; this depends on how close the dressed states of |ei) 
and |e2) are shifted into resonance with the driving (£7) by their cou- 
pling (5). 



above, we obtain a single effective Lindblad operator 



1 



^ = (V7|0)<1|) T |1)<1| )[-\l)(0\ 



A 







2A - i-f 1 

The effective scattering rate is thus given by 



7cff = |(o|X cff |o)| 2 



7 n 2 



4A 2 + 7 2 ' 



(40) 



(41) 



The effective Lindblad operator L c g describes Rayleigh 
scattering, i.e., elastic scattering of incident laser photons by 
the transition |0) <-> Seen from the atom this effect will 
contribute to the effective dynamics of state |0) not as a decay 
but as a dephasing of potential coherent couplings to other 
states. 

Above we have given a rather simple example which could 
also be easily solved without these techniques. For more com- 
plicated situations the formalism developed here is highly use- 
ful. In particular, in the following we consider a more compli- 
cated scheme involving four levels which is relevant for dissi- 
pative state preparation by engineered decay. 



F. Engineered decay 



Here, the goal is to prepare a certain desired steady state as 
the outcome of the evolution of the open system. This is done 
by tailoring its dissipative dynamics. To this end we use the 
effective operators of Eqs. <[3j and (|4]) to take nontrivial inter- 
actions between the excited states into account and to identify 
effective decay processes of the system. By an appropriate 
choice of the system parameters these processes can be engi- 
neered by tailoring the "propagator" -ff^H- 

An example of this is depicted in Fig. [3] showing a four- 
level system consisting of two ground states \gi), \g 2 ) an d two 
excited states lei), Ie2). The excited-state Hamiltonian 



A|ei)(ei| + <5|e 2 )(e 2 | + <?(|e 1 )(e 2 | + |e 2 )(ei|) (42) 



contains detunings of A for \ei) and 6 for |e 2 ), respectively, 
and a coupling with a strength of g between |ei) and |e 2 ). We 
assume that there are no processes between the ground states 
so that H g = 0. The weak classical driving described by 



V 



9. 



(\9i)(ei\ + \ei)( 9 i\) 



(43) 



drives the system between \gi) and |ei) with a Rabi frequency 
of ft. We assume |ei) to decay to \g\) at a rate of 7 and |e 2 ) 
to | g 2 ) at a rate of k, represented by the Lindblad operators 



L 



7 — 

i, = 



(44) 
(45) 



/llffiXeil, 
/ «l5'2)(e 2 |. 

The non-Hermitian Hamiltonian is then given by 

H NH =A| ei )(ei| +5\e 2 )(e 2 \ + g{\e x ){e 2 \+ H.c.) , (46) 

with complex detunings 5 = S — ~ and A = A — ^. Us- 
ing Eqs. ^ and Q we obtain the effective Hamiltonian and 
Lindblad operators: 
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"T 6 \SA- g 2 
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2(<5A - g 2 ) 

T/JigQ 
2{g 2 - S A) 
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(47) 

(48) 
(49) 



We note that the effective Hamiltonian H c s only contains a 
shift of \gi). The effective decay process L£ s effectively pre- 
pares the ground state \g 2 ) from \gi) at a rate of 



«cff = \(g 2 \L* s \ gi ) 



ng 2 n 
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4| ff 2 - 5A| 5 



(50) 



The other process L^ ff is a dephasing of \gi) with a rate 



lcS ^ \( gi \L2 s \ 9l )f 



jisfn 2 

A\SA - q- 



(51) 



Our formalism allows us to engineer schemes that serve the 
purpose of dissipative state preparation, such as Refs. EE). 



The strength of the effective Lindblad operator concept is ob- 
vious: We immediately derive the effective pumping rates and 
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FIG. 4: (Color online) Effective and full time evolution of a system 
with engineered decay. The curves illustrate the preparation of the 
ground state | (72) (blue lines, starting from 0) by decay of the ground 
state | gi) (green lines, starting from 1). For weak driving < 7/5 
the effective operators (thick solid) match the full dynamics (dash, 
fl — 7/ 10) very accurately. With increasing driving Q the assump- 
tion of perturbative excitation is no longer valid and the effective 
evolution deviates from the full dynamics (dash-dot, = 7/2; dot, 
= 7). For the simulations the system parameters 7 = k = g/10 
and the optimized detunings 5 — g 2 / A and A — g \f^]~k were used. 



dynamics of the ground states from the initial operators. If one 
desires to optimize the preparation of \g 2 ) from \gi) which 
happens at a rate of k c s, this can be realized by an appropri- 
ate choice of the system parameters, A, 6, and eventually, g. 
Let us assume that 7, k, and g are fixed, that the coupling 
g is strong, g ^ 7, n, and that the detunings A and 5 are 
adjustable. Then the optimum is reached by adjusting the de- 
tunings to 5 op t — g 2 1 'A and A opt = gy/j/n, which leads to 
a maximized effective decay rate of 



„°pt 

v eff 



o 2 
87' 



(52) 



To compare the effective with the full dynamics, we perform 
simulations of the evolution of the system by numerically inte- 
grating the master equations ([T]) and Q. The resulting curves 
are plotted in Fig. [4] Here we show dissipative prepara- 
tion of I (72) from |<7i) for the optimal choice of A and 5 as 
given above. We note that for weak driving f2 (solid lines, 
SI = 7/IO), the curves of the effective and full dynamics ex- 
hibit excellent agreement. For stronger driving fl (dash-dot, 
ft = 7/2; dot, O = 7), the effective dynamics exhibit in- 
creasing deviations. These result from the breakdown of the 
assumption of weak driving used to derive the effective oper- 
ators. 



1. Interpretation and application of the inverse non-Hermitian 
Hamiltonian iJ^H 

In general, a good physical understanding of the effective 
decay mechanisms of an open quantum system is desirable. 
Even more so, it is essential for developing dissipative state 
preparation schemes such as the ones in Refs. Here, 
engineered decay processes are tailored to prepare a desired 
steady state. In the following, we discuss the physical mean- 
ing of the elements of the inverse non-Hermitian Hamiltonian 
of the excited-state subspace -ff^H f° r trie exam P ie at hand. 
We find that ff N ^ can be written as 



off l e l 



)(ei\ + ~5j\e 2 )(e 2 \ + 



+ g e s(\ei)(e 2 \ + \e 2 )(e 1 \ 
having defined the quantities 



»cff 



= A — g 2 /~S 
~S cS = l/(e 2 \H^\e 2 )=S-g 2 /A 
g cS = l/(ex\H^\e 2 ) = g~6 ■ A/g. 



(53) 

(54) 
(55) 
(56) 



Each of the latter quantities can be seen as an effective com- 
plex detuning (A e g, d e g) or coupling (g e s) of the excited 
states. contains their inverses which act as "propaga- 

tors" for the effective operators of Eqs. Q and Q and there- 
fore govern the strength of the effective processes. Thus, we 
can also express the effective decay rates of Eqs. pO} and 
( |5T| ) in terms of the effective complex energies and couplings 
of the excited states: 



Keff 



7eff 



K n 2 

4|A cff |2 



(57) 
(58) 



We note that our above choice of 6 = S opt minimizes |A e ff | 2 , 
|<5 e ff| 2 , and |g c ff| 2 - Physically, this corresponds to the case 
where the driving V is in resonance with the lower dressed 
state of the excited states |ei) and |e 2 ) (or the upper dressed 

state for 5 = — respectively), as can be seen from Fig. 3 
Accordingly, the absolute values of the propagators of 
in Eq. ( |53~| are maximized under this choice, resulting in an 
enhanced decay from \gi) to \g 2 ). A more detailed account on 
engineered decay mediated by dressed excited states is pro- 
vided in Ref. Q. 



IV. EXTENSIONS OF THE EFFECTIVE OPERATOR 
FORMALISM 

A. Nonperturbative ground-state coupling 

The formalism of Eqs. |2j) - (|4]) was derived assuming that 
the ground-state couplings H g are much weaker than those of 
the excited states contained by J?nh- Under this assumption 



8 



of a perturbative ground-state coupling it was possible to ne- 
glect the effect of H g on the effective processes. 

For strong interactions between the ground states, the action 
of H g can no longer be ignored in the effective processes so 
that the accuracy of the effective dynamics of Eqs. Q and (|4| 
will decrease when H g approaches 7?nh- We will now show 
how to overcome this drawback by diagonalizing the ground- 
state Hamiltonian H g and including its action in the effective 
operators. 

We can build on the dynamics of the separate subspaces as 



given by Eqs. (25 1 and (27 1, derived without any assumption 
about the strength of H g . In contrast to the above derivation, 
we can no longer assume 0(t)P g ~ lP g . Since H g and Hnh 
do not couple the ground and excited subspaces, we can sep- 
arate the evolution operator into one part for each of the sub- 
spaces: 



0(t) = P g O g {t)P g + P c O c (t)P c 

— g—iHgt p _|_ e —iH NH t p 



We assume that H g can be diagonalized 



(59) 



(60) 



with dressed state energies Ei and a projector Pi = for 
each ground state I. Accordingly, we decompose the pertur- 
bative excitations V+ with respect to the ground states, 



(61) 



where we have defined Vi 



V+Pi as the excitation from 



ground state I. Given the different energies of the dressed 
ground states, the effective evolution is now no longer identi- 
cal for ground states with different energy Ei . This additional 
complication in the integral is taken into account by introduc- 
ing the sum in Eq. ( |61) , For instance, for the first term of Eq. 
p7| > we have 

h = P g V(t) V /' dt' fyt)p, (< V)fl (62) 

, JO 



With Eqs. (59 1-(61 1 this term becomes 



h =p g 6 g 1 {t)V-d c {t) V / dt' e t( - H ™- El)t 'vlp (0 \t) 

i J o 

« iP s 6-\t)v. ]T (hnh - EtY 1 vld e (t)pW(t). 



The integration is carried out similarly for the other terms in 



Eq. (27 1 (the "sandwich" terms can be neglected for the same 
reasons as above). Transforming back into the Schrodinger 
picture and arranging the terms as in Eqs. ([3]l and Q we find 
for the effective operators, including nonperturbative ground- 




FIG. 5: The three-level Raman system. Two ground states |0) and 
1 1) are driven up to an excited state je) with different detunings Ao, 
Ai and Rabi frequencies Qo, Hi- |e) decays to the ground states via 
spontaneous emission at rates of 70 and 71 . Effects originating from 
nonperturbative interactions between the ground states (indicated by 
f22) can also be taken into account by our extended formalism. 



state evolution, 



V 



E 



H 



(i) 



NH 



Vi + B.c. 



r(0 
'NH 



vi 



H s , (63) 
(64) 



We see that in order to apply this formalism rather than Eqs. 
( 3]4 1 we replace the general inverse non-Hermitian Hamilto- 
nian i?NH by initial-state dependent propagators 

1 = (ff NH - £j) 1 (65) 

for each ground state I. In doing so we obtain the accurate 
effective dynamics in the presence of nonperturbative ground- 
state coupling. In Ref. [7| we have used this technique to 
investigate rapid preparation of entanglement by engineered 
decay. Below we consider a simpler example of a three-level 
Raman system. 



B. Example: The three-level Raman system 

A three-level system in Raman configuration is a widely 
used quantum system so that the understanding of its effec- 
tive processes is highly relevant. In particular, Ref. [9| deals 
with its effective dynamics in the absence of decoherence in 
great detail. In the following we give a description of the ef- 
fective dynamics of a three-level Raman system that includes 
dissipation. 

As illustrated in Fig. [5] the system consists of two ground 
states |0) and |1) and an excited state |e). Coherent driving 
of the transitions |0) -R> |e) and |1) f-> |e) is facilitated by 
two fields, generally with different detunings, A and A 1; and 
strengths, SIq an d in a time-independent frame the system 
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is described by 



-A 



AillXll, H e = 



V= ^(|0><e| + |e)<0|) + ^(|l)<e| 



(66) 

|e)<l|), (67) 



assuming an arbitrarily strong nonperturbative H g . The time- 
independent formulation of Eq. ( 66 1 with detunings Ao and 
Ai assigned to the ground states allows us to use the formal- 
ism of Sec. IV A Alternatively, we could load the time de- 
pendence on the fields and use the formalism to be presented 
in Sec. |IVC| below. 

and 



ground states 
the Lindblad operators 



Decay from the excited level |e) into 
1) at rates of 70 and 71 is described by 



^7,0 



^7,1 



V^|0)(e| 
7l|l)(e|. 



(68) 
(69) 



The non-Hermitian Hamiltonian can be divided into two parts, 
denoted by the initial state of the exciting field 



(o) 

NH 



H 



(i) 

NH 



A 



Ai 



' y 

i'y 
~2 



|e)<e| = A |e)(e 
e )(e|^A 1 |e)(e| 



(70) 



(71) 



In the last step we have once again assigned complex energies 
that combine the real detuning and the imaginary decay of the 
levels. 

Using the non-Hermitian Hamiltonians of Eqs. (70 1 and 
together with Eqs. (63i-(-67i we obtain the effective 



(71 



Hamiltonian 



H eS = A 



0"o 



A1- 



4A 2 + 7 2 

Ai^ 

4Aj + 7 2 



(72) 



liXil 



(A + Ai)ft fii 



-|0)(l|+ff.c. 



v 8(Ao-i7/2)(A 1 +i 7 /2) 1 

where we have defined an overall decay rate of 7 = 70 + 71- 
We note that despite the complex terms the effective Hamilto- 
nian is Hermitian. Besides two shift terms similar to the one 



in Eq. ( 39 1, it contains an effective two-photon transition be- 
tween the two ground states with an effective Rabi frequency 

_ f o _ (Aq+Ai)O Oi 

Ol Si cff _ 8 ( Ao _ i7 /2)( Al +i7/2)- 

In the absence of dissipative processes (7 = 0), this ef- 
fective Hamiltonian equals the time-averaged ground-state 
Hamiltonian of Gamel and James [9], viz., 

ffe ff = -^ + ^(|0Ki| + |iK0|)- 



o 2 

s 'o 



(|0)(1| 



4A 4Ai / 

"4a ( ; |Ua0| "4^ 1 

Furthermore, we derive the effective Lindblad operators 

/To^o lriWril , y/lo^i 



f7.0 



L 



2(A - vy/2) 
- .1 = VtT^q , 
off 2(Ao-n/2)' 



|0>(0| 



2(Ai-i 7 /2) 



2(Ai - i 7 /2) 



|0)(1| 



(73) 



(74) 



(75) 



a) 1.0 



c 
o 



3 
Q. 
O 
Q. 



0.9 



0.8 



0.2 



0.1 



0.0 



A A 

:V\ 


/ 


VA 






— fP ) (P ) 

» 1 'full v 1'eff 

(P ) (P ) 

v 0'full v 0'eff 


" A A\ 






\ \ 
\ \ 





0.0 5.0x1 2 1.0x10 3 1.5x10 3 2.0x1 3 

t (units of 1/A) 



b) 



c 
o 

Q. 
O 
Q. 




0.0 4.0x1 3 8.0x1 3 1.2x10" 1.6x10 4 2.0x1 4 
t (units Of 1/A) 

FIG. 6: (Color online) Comparison of effective and full evolution of 
a three-level Raman system. Curves obtained from numerical inte- 
gration of the effective master equation (dashed) agree with results 
from the full master equation (solid) both (a) in the absence (7 = 0) 
and (b) in the presence of dissipation (7 7^ 0). The effective opera- 
tors are found to model the slow dynamics of the two ground states 
0) (blue, starting from 0) and |1) (green, starting from 1) very ac- 
curately, averaging out the fast oscillations (inset in a). For the sim- 
ulations we used the parameters fio.i = A/10 (A = Ao + Ai): 
(a) Ai - A = A/1000, 70,1 = 0, (b) Ai - A = A/100, 
70,! = A/10. 



Besides one loop-term for each of the ground states in this 
setup, these operators contain effective decays from either 
ground state to the other, the strength of which is given by 



7^1(11^ 



•7,1 1 



7eV° = 



l^|l)| 2 =, A2 



7i "§ 
4A 2 + 72 ' 

7o«i 



4A 2 + 7 2 • 



(76) 
(77) 



We note that depending on the relative strength of the effec- 
tive quantities f2 c ff and j e g, the resulting effective dynamics 
will either be governed by coherent or decoherent behavior. 
This is visualized in Fig. [6] where we have plotted simulated 
curves of the full and effective dynamics obtained by numer- 
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ical integration of the master equations ([T]i-(j2]i. In Fig. |6|a) 
we show a purely unitary case, and in[6jb) a mixed case with 
both coherent and dissipative processes present. 

In the purely unitary case (7 = 0) shown in Fig. |6ja) we 
see that the populations of the two ground levels exhibit Rabi 
oscillations at a high and a low frequency. The high-frequency 
oscillations correspond to the virtual excitation of the excited 
state |e). These oscillations are explicitly excluded from the 
formalism developed here and are therefore not present in the 
evolution with the effective operators, as can be seen from the 
inset in Fig. |6|a). Nevertheless, the formalism captures the 
slow dynamics of the ground states. 

A case including dissipation (7 7^ 0) is shown in Fig. |6|b). 
Here we see that even for Rabi oscillations sweeping almost 
the entire population between the ground states, the effective 
dynamics match the time evolution of the full master equation 
with very high precision. For large times t the oscillations are 
damped out and the dynamics converge into a steady state. 

Finally, we comment on the situation where the ground 
states are coupled by another field of strength f^, as illus- 
trated in Fig. [5] In case this additional interaction is pertur- 
bative it can simply be included in H g , and hence, in H e g, 
without affecting the other terms. In the nonperturbative case 
the effect of 0,2 on the effective processes is caught by chang- 
ing into a frame in which the ground-state Hamiltonian H g 
is diagonal. From there the formalism of Sec. IV A can be 
applied. 



C. Several perturbations or fields 

In the following we present an extension of our effective op- 
erator formalism to several perturbations or fields Vf, where 
/ denotes the particular field and ojf its frequency. Then we 
can write the perturbations as 



v + e 



B.C. (78) 



The formalism we develop can also be used to include non- 
perturbative ground-state coupling as in Sec. |IV A| Still, at 
the first glance the assumption of several fields seems prob- 
lematic: so far we assumed a rotating frame of reference in 
which V is time independent. However, our formalism can 
be derived without this claim, starting from Eqs. (25 1 and 



(27 1. As opposed to the previous case, where we chose to 



work in the interaction picture, we now keep the time depen- 
dence in the perturbations. For simplicity we choose a time- 
independent frame with respect to the interactions inside the 

subspaces i^f- = ^f- — 0). Again, for perturbative ground- 
state coupling H g the ground-state evolution becomes negli- 
gible, O g (t) ~ 1. The perturbative evolution of the ground 



states turns into a sum for the different fields /: 

h = P s V(t) [ dt' V(t')pW {t ') Pg = (79) 
Jo 

« P g V-(t)O c (t) I dt' e^^f^' V f + P k p {0 \t) 



« iP s V-(t) £ [Hnu ~ u>f) vl(t)P sP W(t). 
f 

We then obtain the effective operators 
1 



H. 



eff 



V- 



■(*)E(^2) 1 *+(*)+ 



/ 



(80) 
(81) 



with one propagator (-H^h) 1 = (^nh ~ u f) 1 for each 
field/. 



D. General formalism 



In Sec. IV A| we derived an extension of the effective op- 
erators in Eqs. ([3]l and p} for nonperturbative ground-state 
coupling and in Sec. |IV C we provided the extension for sev- 
eral fields. In a situation with both aspects present these two 
extensions can be directly combined. We find the effective 
operators 



Hett — —-, 



'nh 



V_ 



(f,i) 



(t) + H.c. 



f,l 



u 

with field- and state-dependent propagators 



(82) 
(83) 



(84) 



This quantity contains both the information about the initial 
state and the exciting field; thus, it is the most general propa- 
gator expression presented in this work. 



V. COMPARISON TO OTHER METHODS 

In the following we compare the results obtained here with 
other methods from the literature. The formalism we have 
presented here is equivalent to the standard approach of adi- 
abatic elimination in quantum optics (see, e.g., Ref. [2|) and 
is essentially a formalization of it. In most approaches adia- 
batic elimination is done at the level of equations of motion. 
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This procedure can therefore be rather tedious, as it requires 
the derivation of the equations of motion followed by vari- 
ous manipulations of the equations, which are then often used 
to extract effective operators. For comparison our formalism 
works directly at the operators and immediately gives the ef- 
fective operators without reverting to the equations of motion. 

The effective operators of Eqs. Q and Q are the most 
compact formalism presented in this work. For the assump- 
tion of perturbative ground-state coupling and a single pertur- 
bative exciting field these operators match the time evolution 



very precisely. The extended operators of Eqs. ( 63 1 and ( 64 1 



( 80 1 and ( 8 1 1, and ( 82 1 and ( 83 i allow for the same precision 



but also in the presence of nonperturbative ground-state inter- 
actions and several exciting fields. 

As could be seen from the example of the three-level Ra- 
man scheme in Sec. |IV B[ in the absence of decoherence the 
effective Hamiltonian method of James and co-workers EE) 
can lead to similar results as our formalism. However, as op- 
posed to Ref. (9), we do not find any additional decoherence 
terms emerging from averaging over the fast coherent evolu- 
tion of the fields. 

Also, the Feshbach projection-operator method of Refs. 
ifTTl [T9l allows for descriptions of open system dynamics by 
means of an effective non-Hermitian Hamiltonian. We see 
that if we ignore the feeding term J2k ^eSP (Arff)^ m Eq. 



(33 i the evolution is described by an effective non-Hermitian 



Hamiltonian 



it.. 



oil' 



(85) 
(86) 



This Hamiltonian is equivalent to the one of Refs. lfT71[T9l . In 
the language of the Monte Carlo wave function method fl2l . 
-Heff.NH accounts for the "no-jump" evolution of the ground 
states. In contrast to this method, our effective formalism goes 
beyond including effects of non-Hermitian time evolution in 



the effective Hamiltonian, as we also include the feeding term. 
As a result we separate the non-Hermitian effective evolution 
into a Hermitian part with a (Hermitian) Hamiltonian H c r and 
a non-Hermitian part with effective Lindblad operators L\ s . 



VI. CONCLUSION AND OUTLOOK 

We have presented an effective operator method covering 
both coherent Hamiltonian and dissipative Lindblad interac- 
tions. Our effective operator formalism allows us to reduce 
the complexity of an open quantum system considerably by 
restricting its time evolution to an effective master equation 
describing ground-state to ground-state processes. 

Our effective operator formalism is useful for understand- 
ing the quantum dynamics of complex open systems by iden- 
tifying their effective dissipative interactions and by reducing 
high-dimensional evolution to an effective master equation of 
the ground states. More specifically, our effective operators 
can be applied to identify and tailor effective decay processes 
involving coherent driving and naturally occurring sources of 
noise in open quantum systems. In particular, we have found 
the presented operators to allow for the development of phys- 
ical schemes for dissipative quantum computing and dissipa- 
tive state engineering 0171. 
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